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Abstract 

We study the decays of a standard model-like MSSM Higgs boson to pairs of 
neutralinos, each of which subsequently decays promptly to a photon and a gravitino. 
Such decays can arise in super symmetric scenarios where super symmetry breaking is 
mediated to us by gauge interactions with a relatively light gauge messenger sector 
{Mmess ^ 100 TeV). This process gives rise to a collider signal consisting of a pair of 
photons and missing energy. In the present work we investigate the bounds on this 
scenario within the minimal super symmetric standard model from existing collider 
data. We also study the prospects for discovering the Higgs boson through this decay 
mode with upcoming data from the Tevatron and the LHC. 



1 Introduction 



Supersymmetry (SUSY) is an attractive mechanism to stabilize the scale of electroweak 
symmetry breaking against large quantum corrections from unknown high-energy physics. 
However, to be phenomenologically viable SUSY must be broken at low energies. One class of 
models that can achieve this in an acceptable manner are theories of gauge-mediated super- 
symmetry breaking (GMSB) where supersymmetry breaking occurs in a hidden sector, and 
is communicated to the visible sector containing the supermultiplets of standard model (SM) 
particles through gauge interactions with a set of gauge messenger particles [H [21 [3] . 

Gauge mediation leads to experimental signatures that are very different from other forms 
of SUSY breaking, such as gravity and anomaly mediation, when the gauge messengers are 
very light relative to the Planck scale. In this case, the lightest superpartner particle (LSP) 
is generally the gravitino, and is stable in the presence of i?-parity (which we assume in the 
present work)j3 The lightest standard model superpartner will then decay to a gravitino 
and one or more SM states. For example, if the lightest SM superpartner is a mostly-Bino 
neutralino, it can decay to a photon and gravitino. Such decays can be prompt on collider 
time scales for gauge messenger masses below about 100 TeV. The distinctive photon- 
rich collider signatures of this scenario have been studied in a number of works such as 

Refs. HElEllTllHlElllollIIinailS]. 

These previous studies of low-scale gauge mediation focused mainly on a minimal class of 
GMSB models. However, more recent investigations of GMSB scenarios have illustrated that 
a much broader range of superpartner spectra are possible within more general realizations 
of this mechanism [TU |15l [161 [13 UHl HI]. With the Tevatron running and the LHC about to 
start up, it is therefore interesting to consider new and unusual collider signatures that can 
emerge in this context [20] . In the present work we investigate a novel Higgs boson signature 
within the minimal supersymmetric standard model (MSSM) that can potentially arise in 
generalized GMSB models, but is impossible in minimal GMSB scenarios. 

The signature that we discuss arises from the decays of a SM-like Higgs boson to pairs 
of the lightest neutralinos. Each neutralino is assumed to subsequently decay promptly to a 
photon and a stable gravitino, Xi ~^ iG, giving rise to a collider signature consisting of two 
photons and missing energy (Pt)- This decay mode does not occur in standard scenarios 
of gravity-mediated supersymmetry breaking because in that case the lightest neutralino is 
(meta-)stable and can only produce an invisible Higgs signature |2T1 [221 123]. This decay 
channel is also impossible in minimal MSSM GMSB models since the lightest neutralino 
is constrained to be heavier than half the mass of the SM-like Higgs boson. However, 
Higgs boson decays to pairs of unstable neutralinos can potentially arise in generalized 
GMSB constructions. To the best of our knowledge, this decay mode has not been studied 
previously, although related work has treated the decays of the heavier non-SM-like Higgs 
bosons through this channel [23]. Let us also mention that decays of Higgs bosons directly 
to a gravitino and a neutralino have been investigated in Ref. [25], while Higgs decays to 
a heavier and a lighter neutralino were studied in Ref. [26]. Non-standard Higgs decays in 

^i?-parity is one simple way to forbid dangerous operators that can lead to rapid proton decay. 
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generalized GMSB scenarios within singlet extensions of the MSSM have also been discussed 
in Refs. [271 [281 129] . 

The layout of this paper is as follows. In Section [2] we investigate the constraints on this 
scenario from previous searches at LEP as well as direct GMSB searches at the Tevatron. 
We estimate the prospects of discovering this Higgs boson decay mode at the Tevatron in 
Section [31 In Section [H we undertake a similar analysis for the LHC Finally, Section [51 is 
reserved for our conclusions. 



2 LEP and Tevatron Bounds on Light Neutralinos 

We begin by investigating the bounds on light neutralinos that decay promptly to a photon 
and a gravitino from existing particle collider data. Motivated by recent progress in non- 
minimal GMSB models, we consider very general low-energy superpartner spectra (with 
a light gravitino) that do not necessarily fit into the paradigm of minimal GMSB. The 
strongest parameter bounds are found to come from the LEP experiments and the Tevatron 
searches for direct neutralino and chargino production in low-scale GMSB. We also study 
the maximal size of the branching fraction of the SM-like Higgs boson to neutralino pairs 
subject to these constraints. No dark matter constraints are imposed as the dominant dark 
matter component may come from a different sector - examples include axions [301 131] or 
dark matter in the supersymmetry-breaking sector [321 [331 1311 [351 136]. Furthermore, the 
very light gravitino masses considered here need not induce any cosmological difficulties [3Z] • 



The Yukawa and gauge couplings of the lightest MSSM neutralino Xi to the lightest 
CP even Higgs boson and the electroweak gauge bosons are given (in two-component 
notation) by [33 [39] 



where 



y.o^o^o = {cosaN*^ + smaN*,){g'Nl,-gN*,), (2) 

>^.o,o,o ^ ^(liv.p-iivun, 

Here the unitary matrix satisfies the relation A^*A^'^A^^ = diag(m^o, m^o, m^o, m^o) where 

is the neutralino mass matrix in the basis (5°, VT^, if^, if^). Similarly, the matrices U 
and V bi-diagonalize the chargino mass matrix according to U*A4^V'^ relative to the basis 
{W^,H^), and the angle a describes CP-even Higgs boson mass mixing. See [IQ] and 
for further discussion and four-component versions of these formulae. 
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The results of LEP I and II and the Tevatron place bounds primarily on the couplings 
of neutralinos and charginos to the electroweak gauge bosons. These constraints imply that 
in order for the lightest MSSM neutralino to be light enough to allow superpartner Higgs 
decays, it must be mostly Bino since significant Higgsino or Wino components lead to overly 
large gauge boson couplings (and unacceptably light charginos). The Bino, coming from 
f/(l)y, avoids such couplings due to its Abelian nature. This is reflected in the absence 
of A^^ii in the gauge boson couplings listed in Eq. (|2]). However, for the lightest neutralino 
to couple to the Higgs, it must have at least some Higgsino content. From the expressions 
for Higgs-neutralino and neutralino-gauge boson couplings, one sees that two Higgsino (or 
Wino) mixing factors are required for a mostly Bino state to couple to a Z° gauge boson, 
whereas only one is needed to couple to the SM-like Higgs boson. It is this feature that 
will allow the Higgs to decay significantly to a Bino-like neutralino in the MSSM without 
contradicting LEP constraints. 

2.1 Bounds from LEP and the Tevatron 

The most constraining LEP searches for a light neutralino decaying promptly to a photon 
are those looking for di-photon final states. LEP II has searched for di-photons and missing 
energy up to \/i ~ 209 GeV, which limits o"(e+e^ — )■ 77 + ^t) < 10~^ pb [12]. LEP I running 
on the pole also bounds the branching fraction for decays to di-photons and missing 
energy as BR{Z^ — )■ 77 + ^t) ^ 3 x 10~^ |13]. The Z° decay bound was derived assuming 
Z° — )■ uuZ' with Z' — 7- 77 and mz' = 60 GeV. We expect the acceptance for 77 + fix to 
be similar for Z° — )■ XiXi with Xi ~^ lG, and we therefore apply the same bound here. In 
addition to bounds on the neutralino, the lightest chargino must typically be heavier than 
> 103 GeV to be consistent with the LEP data [SI 115] . 

Light neutralinos that decay promptly to a photon and a gravitino are also strongly 
constrained by Tevatron searches for GMSB. These searches performed by the CDF and Z)0 
experiments concentrate on 77 + fix final states. Their results are presented in terms of 
bounds on the specific SPSS GMSB Snowmass point [46], corresponding to minimal gauge 
mediation with tan/3 = 15, N^ess = 1, A = F/M a free parameter, and M^ess = 2A. 
Given their focus on this specific point, their results do not always apply to more general 
scenarios, but we expect their bounds to carry over to the present case. The dominant 
SUSY production mode for the SPSS point at the Tevatron is direct electroweak chargino 
and neutralino creation. In our numerical scans below, we will effectively decouple all 
the superpartners other than the Bino and the Higgsinos, and electroweak production of 
charginos and neutralinos will again be the dominant Tevatron SUSY production mode. 
The primary significant difference between the spectrum we will study and that of SPSS 
is that our lightest neutralino will be much lighter for a given Higgsino mass. This may 
change the decay cascades of the heavier chargino and neutralino states. However, since the 
Tevatron searches concentrate on the 77 + f^T signature, the nature of these cascades should 
not alter the bounds very much. 

The latest CDF GMSB search makes use of 2.6 fb~^ of integrated luminosity [47J. They 
demand well-identified photon pairs, where each photon has > 13 GeV and \r]\ < 1.1. 
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In addition, they apply a cut on the event shape to remove jets mis-identified as photons 
(MetSig), they demand Ht = J2i Pt + > 200 GeV where the sum runs over all hard 
photons, leptons, and jets (suitably defined) to reduce SM backgrounds, and they require 
A0(7i,72) < TT — 0.35. After these cuts, the dominant background comes from electroweak 
gauge bosons producing genuine in the form of neutrinos. In this signal region, 1.2 ± 0.4 
events are expected while none are observed, and we interpret this result as an upper bound 
on the total chargino and neutralino production cross-section of about (Txtot < 20 fb. 

has a similar but slightly different search based on 1.1 fb^^ of data [H]. They require 
well-identified photon pairs such that each photon has > 25 GeV and \r]\ < 1.1. Using 
this di-photon sample, they compare the distribution of observed events to the expected 
background. Using a likelihood analysis based on the signal shape, they find that the total 
chargino plus neutralino production in SPSS must be less than about 80 fb. Note that in 
contrast to the CDF analysis, they do not impose a cut on Hj^. 

2.2 Parameter Scans 

To estimate the maximal size of the Higgs branching fraction to neutralinos in the MSSM, 
subject to the bounds from LEP and the Tevatron on light neutralinos and other superpart- 
ners, we have performed a numerical scan over MSSM parameters using a modification of 
NMSSMtools [49] and PYTHIA6.4 With NMSSMtools we compute the spectrum and 
check the LEP bounds on superpartners and the Higgs, while we use PYTHIA 6.4 to estimate 
Tevatron cross-sections. We have also cross-checked our results obtained in NMSSMtools at 
the qualitative level using CPSuperH [51]. In our scans we fix 

Ms = 700 GeV, M3 = 800 GeV, m^^ = 490 GeV, = 600 GeV, 

mQ^ = muf = 2000 GeV, mzjc = 1970 GeV, A = GeV. (3) 

We also scan over the ranges Mi e [0,80] GeV, fi G [-500,500] GeV, tan/3 G [3,15], 
and m^o G [500, 2000] GeV. Since we restrict our attention to spectra with a light Bino 
and somewhat heavy scalars, these parameters do not obey the relations of minimal gauge 
mediation. However, they may be realized in more general GMSB scenarios where the 
gaugino masses can be taken as independent parameters [H]. In addition, we have adjusted 
the scalar masses to satisfy the sum rules derived in Ref . [16] , although deviations away from 
them can arise in even more general GMSB scenarios^ The large stop and pseudoscalar Higgs 
masses ensure that the lightest Higgs is SM-like and heavier than 114 GeV, even for smaller 
values of tan (3^ We also take somewhat large slepton soft masses to suppress the t-channel 
selectron contribution to LEP neutralino production. We have also studied lower stop and 
Higgs masses, but as we will discuss below, this does not substantially increase the region of 
allowed parameters. 

In the left panel of Fig. [T]we show the values of the branching fraction B{h^ — )■ XiXi) 
obtained in our scan, both before and after imposing bounds from the Tevatron (cr^t^t < 

^ We do not specify here the origin of the /i term, but it is important to keep in mind that mechanisms 
to generate it may induce violations of the sum rules in the third generation or contribute to v4-terms [52j . 
•^Note that the Higgs mass predicted by NMSSMtools has a theoretical uncertainty of a few GeV. 
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20/6). With only the LEP bounds, branching fractions as large as 0.45 are possible for 
m^o > Mzjl- For lighter masses, the constraint from Z° decays limits the Higgsino content 
of the light neutralino, and in turn its coupling to the SM-like Higgs boson. As a result, 
the allowed branching ratio of the Higgs to neutralinos falls rapidly when m^o < 
When the additional bounds from Tevatron GMSB searches are included, the maximal Higgs 
branching fraction to neutralinos is reduced to about 0.15. Branchings of this size are too 
small to significantly affect the LEP-II limits on the SM Higgs mass of about 114 GeV 
derived from searches for Higgs-strahlung off a boson with — > 66 [ 53] . 

In the right panel of Fig. [1] we show the (unboosted) decay lengths cr of the product 
neutralinos as a function of the neutralino mass for a supersymmetry-breaking scale of 
= 50 TeV. We show only points for which BR(h^ — )■ XiXi) > O-l- These decay 
lengths were obtained using the standard result for minimal gauge mediation |54[ [55] 

cr = 487r ^' — ^, (4) 

where |Pi^| = \NiiCw + Ni2Sw\, and ms/2 = \F\/^/3Mpl ~ 0.6eVfl This value of 171^/2 is 
cosmologically safe [2Z] , but still large enough that direct and Higgs decays to gravitinos 
are negligible [251 EE]. Note that the boosts here are typically of order unity because, for 
regions of parameter space where B{h^ XiXi) > O-l; ""^x" '^^^ significantly lighter 

than m/j/2 due to the bound from Z^ decays. Among the allowed points shown in Fig. [H we 
find that they all have decay lengths less than 2 cm, about the limit of what can detect 
using electromagnetic calorimeter (ECAL) pointing. The resolution (based on timing) at 
CDF is somewhat worse. In both cases the small decay length before the photon is emitted 
is not large enough to interfere with the photon reconstruction algorithms. We concentrate 
on such "prompt" photons in the present work, but it would also be interesting to look at 
the case of displaced photons. 

To see the origin of the reduction in the neutralino branching fraction when the Tevatron 
bounds are included, we show in Fig. |2] the values of fj, and tan 13 for the LEP-allowed scan 
points with BR{h^ — )■ XiXi) > 0.1 as well as the range of Tevatron production rates. The 
constraint imposed on the Higgs branching ratio restricts the set of points appearing in the 
^-tan/3 plane. Tevatron bounds require larger /i and tan/3. However, lower values of yU and 
smaller values of tan /3 are needed to obtain a significant Higgs-neutralinos coupling. Let 
us also point out that our scans include both positive and negative values of yU, but only 
positive values generate a significant Higgs branching fraction, greater than 0.1. This arises 
because, having fixed Mi > 0, positive /i leads to more Bino-Higgsino mixing. From Fig. [2] 
we also see that the Tevatron direct production cross-sections of neutralinos and charginos 
become unacceptably large for /i < 250 GeV. This is simply the result of the charginos and 
heavier neutralinos derived from the Higgsinos becoming light enough to be produced at the 
Tevatron without too much kinematic suppression. 



The precise values of cr and TO3/2 will differ in more general GMSB scenarios, but we expect their 
values to be qualitatively similar in many cases. 
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Figure 1: Higgs boson branching fraction into neutralinos — )■ 



(left) and neutralino 



decay length to photon plus gravitino (right) as functions of the neutralino mass. The green 
points in both plots are consistent with LEP bounds while the red points are also consistent 
with the Tevatron search bounds for neutralinos that decay promptly to photons. All points 
in the plot on the right have BR{hP — > XiXi) > 0.1. 



3 Di-photon Searches at the Tevatron 



To investigate the prospects for discovering a SM-like Higgs boson at the Tevatron through 
neutralino decays leading to prompt photons, we have generated parton-level (but including 
initial- and final-state radiation) events in PYTHIA 6.4 ^50j for a particular MSSM sample 
point. The parameter values for this point are 

Ml = 50 GeV, /i = 300 GeV, tan/3 = 5.5, m^o = 1000 GeV, (5) 

with all other soft parameter taken as in our previous parameter scans. For these values 
we find BR{h^ XiXi) - 0.11, mh ~ 114.7 GeV, m^o ~ 46.6 GeV, as well as a total 
leading-order chargino/neutralino Tevatron production cross-section of cr^ ~ 7.2/6. This 
point appears to satisfy all current experimental bounds (not including the DM density). 

Very significant SM backgrounds to the signals we are interested in come from jets and 
electrons mis-identified as photons. We do not attempt to model this detector-dependent 
background. Instead, we simulate signal events while applying cuts and scaling by efficiency 
factors to allow for direct comparisons to the backgrounds tabulated in the most recent 
published searches for — )■ 77 [57] and GMSB [48j. Both searches select events 
containing a pair of photons, each with > 25 GeV and \ri\ < 1.1. In addition, the 
candidate photons must be isolated and pass a set of photon identification requirements 
which differ slightly between the two analyses. The GMSB analysis also considers cuts 
on missing transverse energy (.^r) ^8] . 

The scaling factors we apply consist of a K-factor to account for higher-order corrections 
to the cross-section, as well as a detection efficiency factor e. For Higgs production, we choose 
the K-factors such that the cross-sections in PYTHIA agree with published N^LO results. 
We consider production through gluon fusion (GF) [58l[59], VT/Z-associated modes [60| [61]. 
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Figure 2: Allowed MSSM parameter regions in the /i— tan/3 plane (left) and the corresponding 
Higgsino fraction (|A^i3p + \Nu\^y/^ of Xi (right) for which BR{h^ -> XiXi) > 0.1. The 
green points satisfy the LEP bounds, while the red points also avoid the constraints from 
GMSB searches at the Tevatron. 



and vector boson fusion (VBF) [62], and we obtain scaling factors of Kgf '■ K^/z : Kvbf = 
4.7 : 1.2 : 1.0. The K-factor we use for chargino/neutralino production at the Tevatron is 

The efficiency factors e quoted in Refs. [571 SH] include both the acceptance probability 
A that the signal passes the given cuts, as well as the reduced efficiency e to reconstruct the 
fraction of signal events passing the cuts. We assume that these factorize according to 

6 = Ae. (6) 

While the acceptance A clearly depends on the process under study, we assume further that 
the reduced efficiency e is constant for all events that pass the cuts. To extract the reduced 
efficiencies e, we estimate the signal acceptances for di-photons from a 120 GeV Higgs boson 
or from the SPSS GMSB Snowmass point (A = 100 TeV) by simulating these signals and 
imposing the relevant cuts within PYTHIA. We find Ah — 0.42 and Agmsb — 0.53. Given 
the values Eh — 0.20 for di-photons from a 120 GeV Higgs [57j, and Sgmsb — 0.20 for 
the SPSS point [4SJ, we deduce that Sh — 0.47 and Sgmsb — 0.3S. Presumably the difference 
between these values is due to the shghtly different photon identification requirement used 
in the two analyses. 

In Fig. Owe show the di-photon invariant mass distribution of Tevatron events originating 
from the Higgs, both with and without a cut on missing energy. We also show in this figure 
the di-photon mass distribution from the endpoints of chargino and neutralino cascades. 
All events are required to have at least two photons with > 25 GeV and < 1.1. 
The distributions for which a f^T cut was not imposed have been scaled to be directly 
comparable to the backgrounds tabulated in the di-photon Higgs search in Ref. [57], 
while the distributions with > 30 GeV are scaled to be comparable with the GMSB 



S 




Figure 3: Di-photon signal events at the Tevatron from the Higgs boson and 
chargino/neutralino cascade decays. The distributions without a f^T cut are scaled by K- 
factors and efficiencies to allow for direct comparison with the backgrounds of Ref. [57], while 
the distributions with the imposed cut f^T > 30 GeV are scaled to make them comparable 
to the backgrounds in Ref. [18] . 

search [IS] El No detector smearing effects have been included in this plot. 

Fig. 13] illustrates that di-photons from the Higgs dominate strongly over those from 
chargino and neutralino cascades in the low m^^ region. Di-photons from Higgs decays to 
neutralinos are very spread out relative to direct di-photon decays, which appear as the 
sharp spike around m^^ ~ 115 GeV. The broad distribution of di-photon events from Higgs 
decays to neutralinos is also seen to fall off sharply near 90 GeV. This arises from a kinematic 
endpoint in the invariant mass distribution, 

m^f < , (7) 

rrih - J ml - Arril^o 

with equality occurring when both photons are emitted parallel to the direction of the outgo- 
ing neutralinos. In principle, a measurement of this endpoint together with an independent 
determination of the Higgs boson mass from the location of the diphoton resonance would 
yield the neutralino mass. 

^Note that the QCD background in the Higgs search is estimated from a sideband analysis that cuts out 
the window m;i±15 GeV. Di-photons from — XiXi '^iU ^-l^o contribute in this sideband region. However, 
the signal is much smaller than the background within any one bin, so this effect will be tiny. 
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Comparing our simulation to the background estimate in Ref . [57] , the signal from — )■ 
XiXi appears to be too small to be observable. From a pure counting perspective, there 
are about S = 7 signal events with 50 GeV < m^^ < 90 GeV relative to about B = 8000 
expected for background (3 and 2000 for 70 GeV < m^^ < 90 GeV), yielding a significance 
S/^/B < 0.08 (0.07) for 1 fb-^ of dataj^ Increasing the integrated luminosity to 10 /6 is 
still not enough to yield a visible signal. Even with much more luminosity, the very small 
value of S/B in this channel will likely make it impossible to observe. 

To reduce the SM background, a simple option is to impose a cut on missing energy ^t, 
as considered in the searches for gauge mediation with prompt photons [H]. Adding a cut 
of > 30 GeV to the signal detection requirements described above, we find a reduction 
in the — )• XiXi ~^ 'j'yGG signal by more than a factor of 2. In comparison, the SM 
background is reduced by a much larger factor, as can be seen in Fig. 1 of Ref. [48j, which is 
based on an integrated luminosity of 1.1 fb~^. For a fiducial integrated luminosity of 1 fb~^, 
we find 2.7 signal events with m^^ < 100 GeV relative to an expected background of 9.8± 1.0 
(see Table I of Ref. [IH]), corresponding to a statistical significance of S/ \/B = 0.86. Scaling 
up to an integrated luminosity of 10 fb~^, this yields a statistical significance close to 3, 
enough for a 95% exclusion and nearly enough for preliminary evidence. This is comparable 
to or better than estimates for the combined Tevatron reach for a SM Higgs boson with mass 
below rrih < 130 GeV, even with improvements in detection efficiency [63] . 

Our simple significance estimate based on counting will be degraded once systematic 
uncertainties and the small number of total events are factored in. On the other hand, 
there is a great deal of shape information we have not used. The di-photon spectrum from 
Higgs decays to neutralinos is peaked towards m^^ ~ 70 GeV, whereas our expectation 
is that the SM backgrounds will fall off quickly with m^^. Therefore modifying the cuts 
within the m^^-pT plane and fitting to a signal shape in this plane using a two-dimensional 
log-likelihood analysis has an excellent chance of improving the signal significance. This 
technique would also reduce the sensitivity to systematics. We illustrate the distribution 
of di-photon events passing the basic photon identification requirements discussed above 
[p^ > 25 GeV, |?7| < 1.1, and photon identification) in Fig. HI We should also mention that 
somewhat relaxing the requirements on may also help improve the signal significance. 
Our results suggest that this channel would benefit from further study by the experimental 
collaborations. Note that the CDF GMSB search technique as presented in Ref. |17] is less 
useful for finding this Higgs mode since their cut of Ht > 200 GeV removes nearly all of 
this signal. 



4 LHC Searches 

For a Higgs boson of mass less than 125 GeV, the most promising search mode at the LHC 
is inclusive production followed by a decay to di-photons. In the mass region 115 < rrih < 
125 GeV, both ATLAS and CMS should be able to discover a SM-like Higgs with about 

^Our background estimates are based on the very rough parametrization ~ 
(10 fb/ GeV)(Y5^^)-3-^ of the background given in Fig. 3 of Ref. [ST]. 



10 




10 20 30 40 50 60 70 80 90 100 
(GeV) 



Figure 4: Distribution of di-photon events from — XiXi the m^^-pj- plane subject to 
the photon requirements discussed in the text (p^ > 25 GeV, |?7| < 1.1, and photon ID). 
The z-axis units are fb/ {10 GeV)^ after applying the rescalings relevant for the GMSB 
search. 

15 fb~^ of data. Here, we attempt to extrapolate these analyses to the broader di-photon 
spectrum expected from Higgs decays to unstable neutralino pairs. 

The ATLAS inclusive — )■ 77 analysis presented in Ref. [M] requires two well-reconstructed 
photons with p'^ > 40 GeV, 25 GeV and < |r/| < 1.37 or 1.52 < |r/| < 2.37. The total 
detection efficiency for di-photons from a 120 GeV Higgs after applying these cuts (along 
with photon ID requirements) is e ~ 0.36 in the absence of pile-up. As in our Tevatron 
analyses, we assume that e = Ae, where A is the probability that the signal passes the 
and T] cuts, and that the reduced efficiency e is effectively constant. Simulating — )■ 77 
events in PYTHIA, we deduce e ~ 0.59, which is consistent with the general photon detection 
efficiencies discussed in Ref. [M]. When simulating signal events from — )• XiXi; 'we consider 
GF, l^/Z-associated production, and VBF, and scale the production cross-sections by factors 
of Kgf = 3.2 [nSl EH], K^Y/z = 1-3 [Sni EI], and Kvbf = 1 [SS], where the particular values 
are chosen to rescale the PYTHIA values to higher-order estimates. 

In Fig. Owe show the expected Higgs boson signal for the MSSM sample point described 
in the previous section after imposing the cuts and scalings discussed above. In this figure we 
also show the di-photon signal from SUSY cascade events. The dominant contribution in this 
case comes from gluino production (M3 = 800 GeV), but electroweak neutralino and chargino 
production is also significant. The SUSY signal with m^^ < is seen to be subleading 
relative to the Higgs contribution. To evaluate the significance of the — )■ XiXi di-photon 
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Figure 5: Photon invariant mass distributions from Higgs decays and GMSB SUSY events 
with cuts and efficiencies as in the ATLAS di-photon Higgs analysis {p^ > 40, 25 GeV, 
|?7| < 2.37, and photon identification). 

signal, we estimate the SM background by making a power-law extrapolation to the inclusive 
di-photon background in Ref . [M] Q We find that the statistical signal significance S/ y/B 
is optimized if we require 60 GeV < m^^ < 90 GeV, yielding about 224 total signal events 
and 44000 background events, corresponding to S/ \^ 1.1 with 1 fb^^ of data. Scaling up 
by luminosity, a discovery in this inclusive channel based on statistics alone would require 
at least 20 fb~^ of data. This number will be strongly degraded by pile-up and systematics 
given the very low value of S/B, but it should also be improvable with a more clever analysis, 
possibly using a mild pT cut. At this point let us also mention that we have focused on 
ATLAS searches here because the CMS di-photon search has a slightly harder photon px 
requirement of > 40 GeV, 35 GeV [65] which significantly reduces the signal acceptance. 

Both ATLAS and CMS also plan to perform di-photon Higgs searches in the exclusive 
vector boson fusion (VBF) and associated W/h and Z/h channels. These modes have lower 
production cross-sections, but the extra products in the events, a pair of forward tagging jets 
or a lepton, allow for a significant reduction in background. To illustrate the potential power 
of these channels, we focus on the proposed CMS search for di-photons from W/h and Z/h 
associated production [651 ES]. Their strategy is to look for a photon pair together with at 
least one lepton (e or fi). Both photons are required to have p^ > 35, 20 GeV, and \ri\ < 2.5. 

''Our very approxiraate background estimate is — (100/6/ GeV) ( iso'ceV ) based on Fig. 6 of 

the h'^ — > 77 analysis of Ref. [64] . 
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As above, we deduce a reduced efficiency e for the subset of events with photons passing 
these cuts by simulating direct Higgs boson decays to di-photons in PYTHIA. We obtain 
e = 0.56, 0.33 for W/h and Z/h, respectively, along with a ii'-factor of 1.3 for the production 
cross-section to match onto higher-order estimates [60l [61]. Combining these factors, we 
estimate a total of 6.7 signal events per fb~^ of data in the di-photon invariant mass window 
of 20 GeV < m^^ < 90 GeV. The estimated background rates listed in Tables 10.26 and 
10.27 of Ref. [65], without restricting about 28 fb, yielding a statistical significance 

of 5*/ \/B ~ 1.26 with 1 fb~^ of data. This corresponds to a 5 a discovery with about 16 fb~^ 
of data, and the prospects will improve even further once the backgrounds are restricted 
to lie within the signal acceptance window. Note that this is considerably better than the 
estimated significance for the direct di-photon channel [65] due to the larger signal rate. 
More generally, we expect that any search channel that uses relatively weak photon pj, 
requirements along with an additional handle to remove backgrounds will be effective for 
searching for di-photons from — )■ XiXi- Other examples include VBF and tth associated 
production. 

One additional strategy that could prove useful is to make use the ability of ATLAS to 
reconstruct photon "tracks" . Using the ECAL alone, the primary di-photon vertex can be 
identified to within about 1.5 cm along the beam axis [64J. However, a significant fraction of 
hard photons interact with material in the trackers and convert to an e~^e~ pair. A converted 
photon can therefore give rise to charged tracks allowing for a very precise determination 
of the primary photon vertex position. Estimates in Ref. [64j suggest this determination 
could be as good as 5 x 10~^ cm. This is useful in regular Higgs searches for determining 
the location of primary interaction vertex. In the case of neutralino decays to gravitinos and 
photons in low-scale gauge mediation, vertexing the products of double photon conversions 
might allow for a measurement of the neutralino lifetime, even when the neutralino decays 
are otherwise "prompt". 



5 Conclusions 

Decays of a SM-like Higgs boson to a pair of neutralinos, each of which subsequently decays to 
a photon and a gravitino, can arise in generalized gauge-mediated supersymmetric scenarios. 
This channel leads to a signal consisting of a di-photon pair plus missing energy at high energy 
colliders. In the present work we have investigated this possibility within the MSSM, and 
we find that neutralino branching fractions of the Higgs as large as 0.15 can be consistent 
with existing collider bounds. For branching fractions not too far below this upper bound, 
both the Tevatron and the LHC can potentially observe the Higgs boson through this non- 
standard decay channel. 

The most promising Tevatron searches for this Higgs mode look for di-photon events 
with a small amount of missing energy, much like in the existing searches for GMSB. 
Our preliminary analysis suggests that evidence for the Higgs boson can be obtained with 
10 of data. Improvements in detection efficiencies and optimized analysis techniques 
making use of shape information in the m^^ — plane could potentially lead to a discovery. 



13 



although a more careful analysis by the collaborations is needed before a firm conclusion can 
be drawn. 

At the LHC, this unusual Higgs boson decay channel can be probed efficiently in asso- 
ciated production channels that use relatively mild cuts on the photon momenta. Inclusive 
production with a mild cut may also be useful for probing this decay mode, though in this 
case a more careful study of the backgrounds is required. In principle, the neutralino mass 
can be determined from the endpoint of the di-photon mass distribution in this channel once 
the Higgs boson mass is known. Photon conversions may also provide information about the 
lifetime of the unstable neutralino and the mass of the gravitino. 

In the present work we have focused on Higgs decays to neutralinos in the MSSM, but 
even more possibilities can arise in extensions with additional singlets [GTf 168] . In this case 
the Higgs boson can potentially have significant branching ratios to any neutralinos that are 
mostly bino or singlino. Depending on the precise spectrum and pattern of mixing angles, 
various signatures involving photons and missing energy are possible. Another interesting 
possibility occurs if the spectrum also contains a light pseudoscalar oi, allowing for example 
the decays — >■ XiXi with Xi — > aiG. This could lead to collider signatures similar to those 
considered in Ref . [69] , but with modified kinematics and missing energy. We hope to study 
these possibilities in a future work [70j. 

Finally, we would like to emphasize the importance of searching broadly for new physics 
in upcoming collider data. In the present analysis we have found that in many cases cuts 
imposed while optimizing searches for new physics in other channels also cut out the majority 
of the signal in the non-standard Higgs channel we are interested in. Given the wide variety 
of possibilities for new physics beyond the standard model, we feel that it is most prudent 
to cast our experimental nets as widely as practicable. 
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